Event distributions of polymer translocation 
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We present event distributions for the polymer translocation obtained by extensive Langevin dynamics simula- 
tions. Such distributions have not been reported previously and they provide new understanding of the stochastic 
characteristics of the process. We extract at a high length scale resolution distributions of polymer segments that 
continuously traverse through a nanoscale pore. The obtained log-normal distributions together with the char- 
acteristics of polymer translocation suggest that it is describable as a multiplicative stochastic process. In spite 
of its clear out-of -equilibrium nature the forced translocation is surprisingly similar to the unforced case. We 
find forms for the distributions almost unaltered with a common cut-off length. We show that the individual 
short-segment and short-time movements inside the pore give the scaling relations r ~ A'^" and r ~ f~^^ for 
the polymer translocation. 

PACS numbers: 87.15.A-,82.35.Lr,82.37.-j 
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I. INTRODUCTION 

Polymer translocation through a nano-scale pore has been 
under intensive research since the experimental study by 
Kasianowicz et at. ^ and the first theoretical treatment by 
Sung and Park fS] and the related study by Muthukumar jst]. 
The important treatment in Ht] set the valid limits for 
the scaling exponents obtained for the unforced and forced 
translocation. The highly non-equilibrium nature of the pro- 
cess was captured in the first analytical treatment based on 
the tension propagation present in the forced translocation by 
Sakaue f^. Broadly speaking, the present theoretical translo- 
cation research is spanned between these few cornerstones. 

We have previously shown that the forced polymer translo- 
cation takes place out of equilibrium and is mainly determined 
by the balance between the driving pore force and the drag 
force exerted on the polymer on the cis side ivllltl- The drag 
force magnitude was shown to change with the tension prop- 
agation, which shows as an increase in the number of moving 
polymer beads Nm on the cis side. The drag force was shown 
to consequently decrease with Nm as the polymer translocates 
to the trans side. The polymer crowding on the trans side was 
shown to modify this dynamics and be enhanced by increasing 
the pore force. The crowding in turn increases the scaling ex- 
ponent a in T ~ where t is the translocation time and N 
the number of polymer beads, i.e. the polymer length |0,lltl- 
We also showed that a decreases with increasing the pore fric- 
tion igt]. Our description of the process is very closely aligned 
with Sakaue's analytical treatment that has been adopted and 
expanded e.g. in lUOll . Sakaue's treatment and the charac- 
teristics described above have been incorporated in a coarse- 
grained model in ifTTIl . 

The coarse-grained description of the forced polymer 
translocation outlined above says very little about the stochas- 
tic nature of the process. The force balance and tension prop- 
agation fall within the realm of classical mechanics. Although 
they give the correct average characteristics of the forced 
polymer translocation, for example the question of how fre- 
quently the translocating polymer stops and reverses its direc- 
tion and how long segments typically pass through the pore 



"in one go" and how long such transitions take remain unan- 
swered. Questions such as these have relevance for imple- 
menting DNA sequencing. The force balance and tension 
propagation framework can outline the force dependence of 
the observed scaling of the translocation time r with the poly- 
mer length N but cannot describe how this scaling comes 
about in the first place. In this paper our motivation is to char- 
acterize the stochastic nature of the polymer translocation and 
the origin of the scaling by collecting data on the individual 
transitions of polymer segments. 

A fundamental question is how the forced translocation dif- 
fers stochastically from the unforced process. The forced 
polymer translocation is a highly non-equilibrium process 
even for a modest driving pore force. From the pore outward 
on both the cis and the trans side an increasing portion of the 
polymer is lifted out of thermal equilibrium in the course of 
translocation |0,[M[l2tl- The correlations along the polymer 
contour increase due to the polymer being driven out of equi- 
librium. Hence, it is natural to expect that applying the pore 
force may alter the obtained distributions describing the trans- 
fer of polymer segments through the pore. As we will show, 
these distributions also shed light on the way the observed 
scaling relations emerge. We have previously found that the 
waiting-time distribution, i.e. the numbers of transitions n 
that take a time At plotted as a function of At, n{At), obey 
closely the Poisson distribution ||12<1 . The events, i.e. the tran- 
sitions of individual segments were registered with the reso- 
lution of a polymer bond length b. It is clear that the polymer 
translocation for which the translocation time r scales with 
the polymer length cannot be a Poissonian process. One is 
then naturally led to presume that the polymer translocation 
shows additional correlations if observed at smaller length 
scales than b. 

The outline of the paper is as follows: In Section the 
computational model and its relation to experiments is ex- 
plained. In Subsection nil Al we define the events that we reg- 
ister from the simulated polymer translocations. The multi- 
plicative stochastic process is defined and the polymer translo- 
cation as such a process described in Subsection llllBI The re- 
sults are reported and discussed in the remaining parts of Sec- 
tion|III] Finally, in Section|IV]we summarize the main conclu- 
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sions based on our findings. 



II. THE COMPUTATIONAL MODEL 

A. The polymer model 

We use the standard bead-spring chain as a coarse-grained 
model for the polymer. Here, adjacent monomers are con- 
nected with anharmonic springs, described by the finitely ex- 
tensible nonlinear elastic (FENE) potential 
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where r is the length of an effective bond and R = I.Sct is 
the maximum bond length. The shifted Lennard- Jones (LJ) 
potential 



< 2^/^, (2) 



is applied between all beads of distance r apart. The parame- 
ter values were chosen as e = 1.0, a = 1.0 and K — 30/cr^. 
As there is no attractive part in the used LJ potential, the 
model is for a polymer in a good solvent. We apply no bend- 
ing potential, rendering the polymer as a freely-jointed chain 
(FJC), instead of a worm-like chain (WLC). We have previ- 
ously shown that the difference of the two models when sim- 
ulating polymer elasticity in mechanical stretching or in flow 
is insignificant ifisll . 



mer. I is thus the length of the pore. No-slip boundary condi- 
tions applied on the polymer beads hitting the surfaces prevent 
them from entering the wall. The pore is modelled as a cylin- 
drical potential whose center axis is perpendicular to the wall 
and extends through it. The pore diameter is 1.2(7. Inside the 
pore, the cylindrically symmetric damped harmonic potential 
Uh pulls the beads toward the pore axis with force 
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where k = 100, c ~ 1, Vp is the distance of a polymer bead 
from the pore axis, and Vp is its velocity component perpen- 
dicular to the axis. As mentioned in Subsection lllBI the pore 
force f{ri) in Eq. Q is exerted only inside the pore in the case 
of forced translocation. To prevent the pore force from fluc- 
tuating with the changing number (2 or 3) of polymer beads 
inside the pore we calculate the exact fraction of the polymer 
inside the pore and adjust the force magnitude applied on each 
bead inside the pore so that the pore exerts a constant force per 
polymer length. 

There is no explicit adhesion between the polymer and the 
pore in our model. While adhesion is of importance in the 
DNA translocation experiments, see e.g. ITtIi , here we want to 
characterize the basic translocation process prior to embark- 
ing on the effects of relevant additional interactions. The tran- 
sitions and pauses of the polymer segments inside the pore 
hence result from the changes in the force exerted on these 
segments inside the pore. These changes in the force, in 
turn, reflect the conformational changes of the polymer seg- 
ments outside the pore and their changing interactions with 
the Brownian heat baths on the cis and trans sides of the wall. 



B. The dynamics 

The dynamics of the polymer translocation was simulated 
using Ermak's implementation of Brownian dynamics ifTill . 
Accordingly, the time derivative of the momentum of the poly- 
mer bead i is given by 

Mt)^-^P^it)+V^{t)+f{r^), (3) 

where Pi{t), r}i{t), and /(r^) are the friction constant, mo- 
mentum, random force of the bead i, and external driving 
force, respectively. When applied, f{Vi) is constant and ex- 
erted only inside the pore. Velocity Verlet was applied in the 
time integration lllSil . In the present study the parameter val- 
ues in reduced units lfl6ll were as follows: ^ = 0.5 and rii{t) is 
related to ^ by the fluctuation-dissipation theorem. The mass 
of a polymer bead is m = 16. The time step used in the nu- 
merical integration of equations of motion was 5t = 0.001. 
Partly due to the out-of- equilibrium nature of driven translo- 
cation, changing either ^ or m changes the obtained scalingof 
the translocation time with the number of beads, r ~ iV" 

C. The wall and the pore 

The wall comprises two aligned surfaces the distance / = 
36 apart, where = 1 is the Kuhn length for the model poly- 



D. Relation to experiments 

In order to relate the computational force to a physical force 
inside a pore in experiments we need to relate the energy and 
length scale in our model to the physical energy and length 
scale. In our reduced units kT = 1 corresponds to fc^T, 
where ks is the Boltzmann constant, and the physical tem- 
perature T is taken to be 300 K. The correspondence between 
computational and physical length scales can be established 
by taking the polymer bond length b as the Kuhn length for 
the physical polymer In SI units the bond length for our 
FJC model polymer can be obtained as b = 2Ap , where Xp 
is the persistence length, 40 A for a ssDNA lUSll . The pore 
force per bond length in SI units, ftot, is then obtained from 
the dimensionless pore force per bond segment, /, by relat- 
ing fb/ksT = fb/kT. The effective pore force per bond 
of / = 1 thus corresponds to / = 0.52 pN per Kuhn length 
for a ssDNA. Since / = 36, / = 1 corresponds to the total 
pore force ftot = 1-56 pN. If one uses Manning condensa- 
tion leading to drastic charge reduction lfl9[ I20I1 to relate this 
to experiments, in the a-HL pore a typical pore potential of 
w 120 mV would correspond to ftot ~ 5 pN. In the light 
of recent computer simulations the reduction of the effective 
force inside the a-HL pore is not predominantly caused by 
the Manning condensation but by electro-osmotic flow that 
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is driven by the motion of counterions along the surface of 
DNA ll2ll I22I1 . However, experiments indicate that the esti- 
mate based on Manning condensation is in the right order of 
magnitude |[23l . 



III. RESULTS : EVENTS OBSERVED WITH HIGH 
RESOLUTION 

Traditionally, an event in a translocation simulation is de- 
fined as a change of the polymer bead at a reference position, 
here the pore opening on the trans side. In order to discern 
correlations at a finer scale than b we have implemented in 
our model the possibility to register polymer segment motion 
inside the pore with the resolution 6/10. In addition, we rede- 
fine the event as a polymer segment traversed in one direction 
inside the pore. The event-related distributions were extracted 
from at least 1000 individual translocation runs for each pore 
force. The runs ending in the polymers translocating to the 
trans side and the ones where they slided back to the cis side 
were identified. Distributions were determined for all runs and 
the trans and cis cases, separately. 

A. Numbers of events vs traversed segment lengths 

First, we define the event as a polymer segment As that tra- 
verses in either direction inside the pore without pausing. We 
sample these events at constant time intervals tint- Already 
using this somewhat limited definition for an event reveals a 
log-normal distribution for the number of events ue vs As. 
In Fig.[TJa) we show distributions n£;(As) for different pore 
force magnitudes. The distributions for motion toward cis and 
trans are identical for unforced translocation. Increasing the 
pore force / these distributions deviate due to the increasing 
proportion of events toward trans. 

The total time r it takes a polymer to translocate consists 
of forward and backward motion and pauses. In order to find 
how the scaling law r ~ emerges we have to include 
the pauses in the definition of the events and register events 
without the time limit introduced by the time interval tint. 
Accordingly, we define an event to be terminated only when 
the direction of the motion is reversed. The characteristics 
are the same as with the first event definition, only the dis- 
tributions become wider To compare the forms of the dis- 
tributions we normalize them with the maximum number of 
events, fie = ne/rimax- These normalized log-binned dis- 
tributions are shown in Fig. [Ub). For clarity, the separated 
distributions toward cis are shown only for the events defined 
in the first way. 

B. Log-normal distribution and multiplicative stochastic 
processes 

The obtained log-normal distributions suggest that the poly- 
mer translocation is largely dominated by the underlying mul- 
tiplicative stochastic process. For such a process the prob- 



ability of an event Pr composed of the succession of m 
independent events with the probabilities (1 < i < m) is 
given by the product of these independent random variables 
pi") = UtlP^, SO that log pi") = Y.ti logP- With m 
large, logPi™-' becomes a normal and, hence, pi™-* a log- 
normal distribution due to the central limit theorem ll24ll . 

In the case of polymer translocation this multiplicativity 
arises naturally, since we observe traversed polymer segments 
As„. Here, n indexes the registered segments. The traversed 
segments As„ toward trans result from m consecutive tran- 
sitions of distance s; toward trans. An analogous statement 
applies, of course, for transitions toward cis. Denoting the 
probabilities of these individual transitions as pi (1 < i < m), 
we see that the probability of m such consecutive transitions is 

Pr = YiiLi Pi- Hence, the polymer translocation is seen to 
be of the simplest type of a multiplicative process. Since the 
traversed segments As„ = msi, the probability distribution 
P(As„) is of the same form as Pr"^\ Taking the continuum 
limit, the enumerable set {As„} is replaced by the continu- 
ous variable As. In this limit P(As) is log-normal due to the 
central limit theorem. In what follows we drop the subscript n 
and use the symbol As to denote the lengths of the traversed 
segments registered with resolution b/10. 

There are subtleties in applying the central limit theorem 
to multiplicative processes, see luSll . These are related to the 
effect of rare events. In our case these correspond to large m, 
or As. As we will see, the translocation is dominated by short 
As events, so applying the central limit theorem is valid and 
the obtained log-normal distributions have the natural expla- 
nation described above. 

The probabilities pi appear independently but are depen- 
dent on the stage of the translocation process. Transition to- 
ward trans gets more probable toward the end of the translo- 
cation process when the majority of the polymer is already 
translocated. Also the crowding of the polymer segments at 
large pore force changes this probability by different amounts 
at different stages of translocation This differing of pi's at 
different stages of the process is characteristic for multiplica- 
tive stochastic processes, see e.g. ll26ll . Due to the additional 
correlations in fluctuations multiplicative stochastic processes 
are typically highly non-equilibrium processes Ii26n27il . which 
we have shown to be the case for the polymer transloca- 
tion 0,11. 

The observed traversed polymer segments As consist of 
successive individual transitions Si measured with resolution 
6/10. The transitions were seen to be slow enough that the 
registered Si = ±6/10. We obtain distributions of the log- 
normal form 

P(As) - cxp[- ln(As/Aso)V2fT2], (5) 
As 

where Aso is the characteristic scale and a the standard devi- 
ation of the variable In As, see e.g. li28il . To further confirm 
that we are indeed sampling events from a distribution result- 
ing from a multiplicative stochastic process we also registered 
each successive transition at times i and collected distributions 
such that even successive transitions of different lengths in 
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the same direction were identified as separate events. Again, 
clean log-normal distributions resulted. 

The forms of the trans distributions are seen to become 
wider with increasing pore force /. This is due to the in- 
creased proportion of long-segment motion at larger /. In 
accordance, the cis distributions become narrower with in- 
creasing /. However, the distributions vary fairly little with 
/. Especially surprising is the small difference between the 
distributions for unforced and forced translocation in spite of 
the fact that forced translocation already at weak pore force 
was seen to be governed by a non-equilibrium force balance 
condition liiSJ. In what follows, we show results for distri- 
butions obtained using the latter event definition. 



C. Numbers of events vs event times 

Next, we look at the distribution of events in terms of 
event times At. Fig. |2] shows the event-time distributions of 
the segments traversed toward trans and cis for / = 1 and 
G {50, 100, 200}. The number of events toward trans de- 
cays as ue ^ e~^^*. The transfer times of the events to- 
ward cis obey ue ^ e~^^* for all At. B « 0.11 for all /, 
whereas A = 0.055, 0.05, and 0.043 for / = 0.25, 1, and 
2, respectively, when N = 200. So, the proportion of events 
taking a longer time increases with /, which is in keeping 
with Fig. [TJb) where the proportion of long-segment events 
are seen to increase with /. This explicitly shows the weak 
increase of the characteristic correlation length with increas- 



ing /. The exponential decay ue 
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is slightly slower 



for polymers of length N = 100 than of iV = 200, see the 
plots for / = 2 in Fig. |2] The dependence of A on even 
for the maximum pore force / = 2 is weak compared with 
its dependence on / and hence is not expected to have any 
effect on the scaling r ^ N". The maximum event times 
and segment lengths are seen to increase only weakly with /. 
In summary, the form of the distribution of events with event 
time varies slightly with / and A^ for events toward trans but 
stays unaltered for events toward cis. 



D. Distributions for constant segment lengths 

For DNA sequencing the important information is the time 
At it takes a polymer to traverse a distance As through the 
pore. Accordingly, we define tB(As) as the total elapsed time 
for transitions of constant distance As. ts consists of multi- 
ple separate events that take slightly different times At. We 
denote the number of such events n£;(Asconst, At). The total 
time spent on transferring segments of length Asconst is thus 
tE = J2'^E {Asconst, At) ■ Ai(Asconst), whcrc the sum runs 
over all registered At for Asconst- In Fig.|3]we show tE for all 
registered discrete constant values of As. In what follows we 
omit the sub text 'const' from As. The obtained distributions 
are of the form 

P{tE{As)) - -^cxp[-lii(As/Aso)V2o-^]cxp(-As/Asc), 
As 

(6) 
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FIG. I. (Color online) Log-binned distributions of events as- 
sorted by the lengths of polymer segments that traverse in one di- 
rection, trans or cis, without pausing and sampled at time intervals, 
or without reversing for different pore force magnitudes /. The 
length of the polymer for all curves is = 200. (a) The number 
of events ue vs. segment length As. ?i_e is an average of num- 
bers of events over several polymer translocations, (b) The num- 
bers of events normalized with the maximum, fiE = riE/nmax 
to compare the distributions. All distributions follow closely log- 
normal forms. The log-normal functions plotted to guide the eye are 
fund = 1/As * exp(-log(As/2)Vl-2) and func2 = 1/As * 
exp(— log(As/2)^/4. The apparent power-law decay at As > 300 
is an artefact due to the log-binning done on very few points there. 
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FIG. 2. (Color online) The event-time distributions ue normal- 
ized with the maximum number of events. Semi-log scale shows the 
exponential decays. See the text for details. 



where Asc is a finite cut-off. All these distributions can 
be fit by keeping Asc and Asq constant and varying cr in 
the log-normal part. The distribution data show fluctuations 
at small As, so we show cumulative distributions for im- 
proved quality. For a continuous variable x the cumulative 
distribution CP of the original distribution P is defined as 
CP{x) = P{x)dx. So, for the discrete variable As we 

have CP(tB(As)) = Ea,=+oo Naturally, afl 

these distributions can be fit with the cumulative form 



CP{tE{As)) - 1 - erf[0n(As/Aso)72(72 + As/AsJ. 

(7) 

Again, the fits can be made by keeping Asc and Aso constant 
and varying cr, which increases with the pore force /. Un- 
like for the unforced translocation, for the forced transloca- 
tion the broader cumulative log-normal distributions can also 
be fit with a decreasing power law and an exponential cutoff, 
As"'^exp(-As/Asco). 

The translocation times can be read off from the maximum 
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FIG. 3. (Color online) (a) The distribution (distr) of event times At 
with segment lengths As for / = and N = 100. The rest are 
cumulative distributions for / = 0, 0.5 and 1, A'^ = 50, 100, and 
200. The cut-off segment length Asc = 95 and the characteristic 
scale Aso — 6 for all fit functions. The dispersion a = 2, 7,and 
9.5 for / — 0, 0.5, and 1, respectively. Also shown is a fit with 
the function pow-exp ~ As^" '' exp(— As/60) for the cumulative 
distribution for f — 1 and A'^ — 50. (b) The cumulative distributions 
of (a) normalized by the maximum value. Lines ^ (As)"""'' and 
~ (As)"'' are drawn to guide the eye. 

values of the cumulative distributions. Statistics is not as good 
as when r is directly measured from and averaged over in- 
dividual simulations. The exponent values a obtained this 
way are a ^ 2.2 for / = 0, and a w 1.4, 1.4, 1.45 for 
/ = 0.25, 0.5, 1, respectively. These approximative values 
are in keeping with our previous results and reported values 
for a in general. The same values for a are obtained by direct 
measurements of r, so this also verifies that the distributions 
were obtained correctly. The dependence of a on / is mild 
compared to some of our earlier simulations. This may be 
caused by our more refined way of taking care that the pore 
force stays strictly constant, see Subsection lllCI The effect of 
fluctuating pore force will be investigated in more detail in a 
forthcoming publication. 

It is noteworthy that for the unforced and forced transloca- 
tion there is a common constant cutoff length Asco ~ 95 that 
does not increase with the pore force or change with the poly- 
mer length A^. The proportion of traversed longer segments, 
however, increases with /, which is seen as a broadening of 
the log-normal distribution. Again, the form of the distribu- 
tion changes very little with / or N, which can be more clearly 
seen from the normalized distributions in Fig|3jb). It would be 
tempting to think that the scaling r ~ N"" resulted from the 
broad log-normal distribution of event times t£;(As). How- 
ever, the exponent values of the power-law segments ( As)"°'^ 
and (As)"^ in Fig|3lb) do not give the correct scaling for r. 



E. How the scaling r ~ A^" comes about 

To more precisely determine the role of events at differ- 
ent length scales we extract events of three different segment 
lengths As = 1, 10, and 50. Extracting only the transitions 
of these three lengths As we obtain distributions of the total 
time t£;(As) depicting the contribution of different individual 
transition times At. These distributions are shown for pore 
force magnitudes / = and 0.25 in Fig. Ufa). The times 
taken up by events of As = 1 are orders of magnitude larger 
than by events As = 10 or 50 for all / G [0, 2]. The peaks 



FIG. 4. (Color online) (a) Distributions of total time t e for constant 
A^ — 200 for transferred segment lengths As £ {1, 10, 50} and pore 
force values f — and 0.25. (b) The distributions for the constant 
pore force f ~ 1 and different polymer lengths A'^. The bottom 
three curves with maxima 13, 35, and 100 are for As = 10 and 
A^ G {50, 100, 200}. The topmost three curves with maxima 130, 
350, and 1000 are for As = 1 and the same polymer lengths. The 
marked maxima give the correct scaling r ~ N", see the text. 
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FIG. 5. (Color online) The distributions of the numbers n_B(As, At) 
for As — 1, and 10 for pore force f — 1. The indicated maxima 
give a = 1.5. 



in the range As £ [0, 20] most clearly visible for / = are 
due to the contribution of events toward cis. With increasing 
/ this contribution diminishes for As = 10 and 50 but not 
for As = 1. This partly explains the dominant contribution 
of very short transitions to the translocation time. For all /, 
almost all the contribution to r comes from events of small 
As and At. As At increases, ts for different As approaches 
the same value. 

The forms of the distributions tsi^s) for constant As are 
identical for all /, including f ~ 0. The maximum values 
of tE for / = 0.25, 0.5, and 1 are max{tE) 360, 180, 
and 105, respectively, when As = 10, and max(<£;) = 4200, 
2100, and 1050, respectively, when As = 1. These values 
give the scaling r ^ /"^ with (3 = 1 for As = 1 and 
/3 < 1 for As = 10. Fig. Ub) shows that maximum val- 
ues of t£;(As) with As = 1 or 10 give the scaling r ~ N", 
where a w 1.5. Since the distributions maintain their form, 
the numbers of events n e (At) for constant As should give the 
same scaling as the total event times tE. In Fig. [5] the event 
time and number distributions for / = 1 and As G {1, 10} 
are shown. The indicated maxima for the number distributions 
indeed give a ^ 1.5 just as the maxima for tE distributions 
indicated in Fig.|4|b). 
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A more accurate estimate for a can be obtained by sum- 
ming up all the total event times or total event numbers for a 
constant As. These give for As = 1: a = 1.45, 1.45, and 
1.47 for / = 0.25, 0.5, and 1, respectively; for As = 10: 
a = 1.45, 1.44, and 1.44 for / = 0.25, 0.5, and 1 are ob- 
tained, respectively. From the direct measurement of the times 
it takes polymers to translocate we obtain a = 1.45, 1.45, and 
1.47 for / = 0.25, 0.5, and 1, respectively. These are exactly 
the same values as obtained from the distributions for As — 1. 
The statistics for / = 0.25 is rather modest, since for such a 
small pore force a vast majority of the simulated polymers 
slide back to the cis side without translocating. Even so, ex- 
actly same estimates for a are obtained by the direct measure- 
ment of T and by summing up all tsi^s = 1). The number 
of short range and time events determine the polymer translo- 
cation even for strongly driven translocation, and the effect of 
mildly changing forms of distributions is negligible. 

IV. CONCLUSION 

In conclusion, we have characterized the stochastic poly- 
mer translocation process with high length resolution. Our 
extensive simulations show that although the unforced and 
forced translocation processes are fundamentally different, the 
pertinent distributions are almost identical in form. The ob- 
tained log-normal forms of the event distributions suggest that 
both processes can be characterized as multiplicative stochas- 
tic processes. The out-of-equilibrium nature and the fact that 
the polymer segment transition probability through the pore 
varies with the process state, i.e. the position of the poly- 
mer with respect to the pore, also fit the general description of 
multiplicative processes. 

Despite the strong non-equihbrium character of forced 
polymer translocation, the variation of scaling exponents in 
the relations r ^ and r ^ do not follow from e.g. 
the pore force changing the pertinent distributions but merely 
from different numbers of events taking place at a fraction of 
the length scale of the model monomer length. Consequently, 



the noise contribution of the thermal heat bath remains sig- 
nificant in spite of the polymer being increasingly lifted out 
of it during the translocation process. The random fluctua- 
tions due to the heat bath restrict the lengths of the individual 
transitions of the polymer segments inside the pore. The ran- 
dom stochasticity exerted upon the segments inside the pore 
result from both the random force exerted directly on these 
segments and the random force resulting from the interaction 
of the dynamically changing polymer conformations and the 
heat bath outside the pore. Although we have shown that scal- 
ing of short transition events give the scaling of the whole 
polymer translocation, the reason for the scaling of the num- 
bers of short transition events remains unsolved. 

As we have shown previously, the most discernible con- 
formational changes are the straightening of the polymer on 
the cis side and the crowding of the polymer segments on the 
trans side iHS]- These changes can be described in a coarse- 
grained fashion even if the stochasticity is omitted. However, 
the data presented here confirms that in spite of its strong 
out-of-equilibrium character the polymer translocation is 
essentially a driven diffusion process when dealing with pore 
force of realistic magnitude. The reported results also show 
the difficulty of controlling the exact position of the polymer 
inside the pore. Although the average translocation can be 
described in a coarse-grained manner based on force balance, 
it will not be of much avail in keeping a polymer segment 
fixed inside the pore. This enhances the importance of the 
ability to precisely control the pore-polymer interaction. This 
interaction should dominate over the stochastic fluctuations 
contributing to the random movement of the polymer inside 
the pore if nanopores are to be used for DNA sequencing. 
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